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-»«. ■ Abstract 

Qh! The four-loop /3-function of quantum chromodynamics is calculated and agree- 

Q H ' ment is found with the previous result. The anomalous dimensions of the quark-gluon 

vertex, and quark, gluon and ghost fields are given for a general compact simple Lie 

group. 
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1 Introduction 

The renormalization group properties of Quantum Chromodynamics were the reason of 
acceptance of this theory as the theory of strong interactions. The central role played 
by the QCD /3-function, calculated at the one- [T], two- j2], three- [Sj and finally at the 
four- loop jlj level, cannot be overestimated in this respect. 

The calculation of the last known, four-loop, term in the expansion of the /3-function, 
was performed by only one group jlj. The authors evaluated "...of the order of 50.000 
4-loop diagrams". These two facts lead to the conclusion, stressed many times (in e.g. 
jHj) that it is not only justified, but also necessary to independently evaluate this quantity. 
The present work is intended to fulfill this need. 

Apart from the /3-function itself, we present two new results. The first one is the 
complete set of MS anomalous dimensions at the four-loop level in the linear gauge and 
with color structures of a general compact simple Lie group. These give, in a compact 
form, the complete set of four-loop QCD renormalization constants. We notice that the 
case of the SU(N) group has been solved in j3] with the assumption, however, that the 
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four-loop /^-function is correct. Our second new result is of a more technical nature, but 
should simplify future renormalization group calculations at this level of the perturbative 
expansion. To this end, we derived the necessary set of divergent parts of the four-loop 
fully massive tadpole master integrals. A purely numeric result was available from [Sj, but 
could not be used for our purpose. 

This paper is organized as follows. The next section presents the methods used in 
the calculation, as well as some efficiency considerations. Subsequently, our results for 
the anomalous dimensions are listed. Conclusions and remarks close the main part of the 
paper. The expansions of the tadpole master integrals are contained in the Appendix. 

2 Calculation 

The calculation of renormalization group parameters, i.e. anomalous dimensions and beta 
functions allows for vast simplifications in comparison to the actual evaluation of Green 
functions with kinematic invariants in the physical region. Dimensional regularization and 
the MS scheme are particularly well suited for this kind of problems, since they make 
it possible to manipulate dimensionful parameters of the theory. In fact, ever since the 
introduction of the Infrared Rearrangement [7], it is known how to set most of them to 
zero and avoid spurious infrared poles. With the advent of the R* operation [S], infrared 
divergences are even allowed in individual diagrams and only compensated by counterterms 
afterward. 

At the four-loop level two techniques seem to be most promising. One is a global R* 
operation [5], where one would set all of the external momenta to zero and also almost 
all of the masses, keeping just one massive line. The spurious infrared divergences are 
then compensated by adding a global counterterm. The advantage of this approach is that 
the whole problem can be reduced to the calculation of three-loop massless propagators, 
for which there exists a well tested and efficient FORM jH] package, MINCER ^U]. The 
disadvantage is that the construction of the global counterterm is not trivial. In fact, up 
to now it has not been possible for the gluon propagator. 

A second technique consists in setting all of the external momenta to zero, but keeping a 
common non-zero mass for all the internal lines fHIU^]. ^ ne advantage of this approach 
is that one never encounters any infrared divergences. One minor part of the price for 
this convenience is the necessity of a gluon mass counterterm. The more problematic part 
is, of course, the calculation of the divergent parts of four-loop tadpole diagrams that 
occur. One way to do this is to generalize the algorithms of ^2] to the four-loop level. 
We, however, decided to use integration-by-parts identities to reduce all of the integrals to 
a set of master integrals depicted in Fig. ^ The divergent parts of the latter were then 
calculated as described in Appendix lAl 

Instead of developing a dedicated software for the reduction of tadpole integrals 1 we 
used our own implementation of the Laporta algorithm ^3] in the form of the C++ library 

1 Such a software has been developed in 
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Figure 1: Completely massive four-loop tadpole master integrals. Dashed lines mark those 
integrals which are needed to one order higher in the e expansion, i.e. up to finite parts 
for this calculation. 
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Table 1: Distribution between the tadpole prototypes of the 203580 different integrals 
occurring in the calculation of the gluon propagator in linear gauge limited to at most one 
power of £. "Denominator powers" denote in fact the total number of dots on the lines. 

DiaGen/IdSolver [T3]. We found it also a good opportunity to study the efficiency of this 
approach on a large scale problem. 

Since the calculation was performed in the linear gauge, the gluon propagator had the 
form 

iDM = ^(-9r + t^)- (1) 

It is clear that this implies that every power of the gauge parameter will lead to more 
powers of the denominators and irreducible numerators in the integrals. A minimal way to 
test gauge invariance at the end is to keep at most a single power of £, which corresponds to 
a first order expansion of the result around the Feynman gauge. Under this restriction, the 
calculation involved a little less than 210.000 independent integrals. The distribution of 
them among the four-loop tadpole prototypes in the hardest case of the gluon propagator 
is given in Tab. H The remaining cases of the quark-gluon vertex, and quark and ghost 
propagators involved about one-third of these integrals and one power of denominator and 
irreducible numerator less for each of the prototypes. 

It turned out that to reduce all of the integrals to masters, it was sufficient to generate 
integration-by-parts identities with up to six additional powers of the denominators and 
four of the numerators for all of the integrals up to seven lines, and with five additional 
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powers of the denominators and four of the numerators for the eight- and nine-liners. The 
total number of solved integrals was then about 2.000.000, meaning a 10% efficiency. About 
37% percent of the integrals turned out to be finite. These could have been eliminated from 
the very beginning by a careful study of divergences. We convinced ourselves, however, 
that this would not allow for lower powers of denominators and numerators in the reduction 
process, unless some very involved procedure were used. 



3 Results 

Since the /3-function is, up to normalization, the anomalous dimension of the coupling 
constant, it is necessary to perform the complete renormalization of some vertex. To this 
end we chose the quark-gluon interaction, mostly because the calculation of 4] involved 
the ghost and gluon instead. 

The renormalization constant of the quark-gluon vertex will subsequently be denoted 
by Zi, whereas the renormalization constants of the quark, gluon and ghost fields by Z 2 , 
Z ?J and Zg respectively. Even though Z% is not necessary for the present calculation, we 
derived it in order to have the complete set of renormalization constants at the four-loop 
level. 

We will not give the renormalization constants explicitly, but instead we will limit 
ourselves to the anomalous dimensions, which are defined by 

where \x is the 't Hooft unit of mass introduced to keep the renormalized coupling con- 
stant dimensionless. Since we use dimensional regularization and the MS scheme, the 
renormalization constants can be expanded as 

i=i e i=i j=i e 

where a s is connected to the QCD coupling constant g by a s = a s /(47r) = g 2 /(167r 2 ). 
Using the fact that the dependence of Z on \i enters only through a s and £, and that the 
renormalization constants of the gluon field and of the gauge parameter are equal, one 
obtains 

I ,a\ dlogZ ^<91ogZ 

(-6 + /3X— 73 (1 - i) = -7, (4) 

where the /3-function is simply equal to the anomalous dimension of a s , i.e. j3 = 7 Qs . This 
implies that 

7 = a,f- = -E< +1 7 (i) , (5) 

Ua s i=Q 
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where now 



7 (0 = + ( 6 ) 

Eq. HJcan be used in turn to reconstruct the original renormalization constant from a given 
anomalous dimension. 

Since the quark-gluon vertex, as well as the quark and gluon anomalous dimensions 
have already been given up to the three-loop level in the linear gauge in 3J, we will not 
reproduce them here 2 but only give our result for the four-loop anomalous dimensions with 
color structures of a general compact simple Lie group. At this point we stress once more, 
that the results have been obtained in a first order expansion around the Feynman gauge 
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2 During the course of our calculation we found full agreement with [5J- Contrary to the four-loop case, 
our three-loop calculation was performed without expansion in the gauge parameter. 
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The first three terms of the expansion of the anomalous dimension of the ghost field 
in the linear gauge can be found in [3]. Even though the results there are restricted to 
the SU(N) group, the values for a general compact simple Lie group can be derived on the 
basis of the fact that only quadratic Casimir operators occur. 

As before our four-loop result is obtained in a first order expansion around the Feynman 
gauge 
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The notation is similar to the one used in Besides Riemann ( functions, the result 
contains the quadratic Casimir operators of the fundamental and adjoint representations, 
C F and Ca, as well as the normalization of the trace of the fundamental representation 
T F 5 ab = Tr (r a T b ^j , where T a are the representation generators, and the number of fermion 
families n/. The higher order invariants are constructed from the symmetric tensors 
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(and similarly for the adjoint representation) and of the dimensions of the fundamental and 
adjoint representations, N F and Na respectively. Using the specific values for the SU(N) 
group 

„ 1 „ N 2 -l n AT df cd df cd N 4 -6N 2 + 18 
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d abcd d abcd N(N 2 + Q) df cd df cd N 2 (N 2 + 3Q) 2 

~^V^ = 48 ' ~^V^ = 24 ' Na = N ( 13 ) 

we checked that Eqs. I71ITU1 are in perfect agreement with [Sj . 

We can now combine the anomalous dimensions to reach the goal of our calculation, 
i.e. the four-loop /^-function. Since, Z as = Z 2 = ZfZ^Z^ 1 , we have (3 = 2^i — 2^ 2 ~ 73? 
and 
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This result, manifestly gauge invariant, confirms jlj. For completeness, we reproduce the 
lower order values, which we have also calculated 
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-C F C A T F n f + 2C F Tpn f , 



Pi = ^C 2 A -^C A T F n f -AC F T F n h 

Po = ^C A -^T F n f . (15) 



4 Conclusions 

We have completed the four-loop MS renormalization program of an unbroken gauge theory 
with fermions and a general compact simple Lie group in the linear gauge. The fact that 
we performed an expansion in the gauge parameter still allows for gauge invariance tests 
in practical calculations, although special gauges such as the Landau gauge cannot be 
chosen. If such need would occur, one may use the SU(N) results from 5\. The correctness 
of the latter relies, however, on the four-loop /3-function. Therefore, the present work was 
not only indispensable to confirm the value of the four-loop /^-function itself, but also the 
correctness of [3] . This goal has been reached. 

It should be stressed that with the accumulated solved integrals, one could easily obtain 
one more term in the £ expansion of the anomalous dimensions. This is due to the fact 
that the complexity of the integrals occurring in the gluon propagator with at most a single 
power of £ is comparable to the complexity of the integrals for the other functions with 
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at most two powers of £. Having solved the latter, the gauge invariant /3-function enables 
then to recover the third term in the £ expansion of the gluon propagator as well. We did 
not perform such a calculation, since from the purely pragmatical point of view, only the 
complete £ dependence would be an improvement. 

A second comment concerns the independence of our calculation. The reader should 
notice that besides basic software as FORM [§], Fermat JH] etc. and our own programs 
|15j . the only external package that we used was the Color package of the FORM 
distribution. On a diagram per diagram basis, we made sure that the color factors produced 
by this package agree with the SU(N) values obtained with the algorithm of [IE]. Our 
confidence was strengthened by the fact that we agreed with |2} on the final results. 
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A Master integral expansions 



Here we give the e expansions of the master integrals depicted in Fig.Qup to the necessary 
order. The notation is similar to that used in MATAD [19 , i.e. we use the integration 
measure 
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with d = 4 — 2e and the constants which occur in the divergent parts of the integrals are 
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where Cl2(a;) = Im(Li 2 (e* z )) is the Clausen function. Notice, however, that we use the 
Minkowski space metric and "Minkowski type" propagators, i.e. l/{k 2 — m 2 ). 
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PRi = j, + ^ + j2 (io + 2 c 2 ) + ^ (20 + 8(2-^3) 

16 

+35 + 20 C 2 -y C 3 + 6C 4 + 0(e). 

The expansions of the irreducible four-loop integrals have been obtained either from 
the finiteness of the same integral with higher powers of denominators or with the help of 
the method described in [T2]. After taking into account the different normalization and 
translating master integrals with dots to integrals with suitable irreducible numerators, the 
above results for the divergent parts are in agreement with the numerical values obtained 
in 0. 

The integrals PR1-PR4, PR4d and PR7 are needed up to constant parts. It is, however, 
sufficient to express the finite part of PR7 through the finite parts of PR4 and PR4d and 
keep them as symbols, since they always cancel from the divergent part of any four-loop 
vacuum integral. 

The reader should notice that we also needed e expansions of two- and three-loop 
tadpoles. These have been obtained in [213 EI]- The values can be read off from our results 
on the reducible four-loop integrals above. 
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